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Stephen S. Stahara 
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I. INTRODUCTION 

This is the final sununary report under Contract No. NASW- 
3184 for the National Aeronautics and Space Administration. All 
of the important results of the research performed under this 
contract have been reported in the open literature, both in 
scientific journals and as technical papers at scientific 
meetings, with appropriate acknowledgements to NASA support. 

This summary report provides a statement of the problem studied, 
a descriptive summary of the most important results, a reference 
list and copies of all publications resulting from the research, 
and a list of all participating scientific personnel. 

2. STATEMENT OF PROBLEM STUDIED 

The problem toward which the research under this contract 
was directed was the modeling of the three-dimensional global 
interaction of the solar wxnd with terrestrial planet magneto/ 
ionospheres. This was to be accomplished based on a coatinuum 
rathei. than particle viewpoint. The theoretical method employed 
is based on an established single fluid, steady, dissipationless 
magnetohydrodynamic model that is appropriate for the calcula- 
tion of supersonic, super-Al f venic solar wind flow past magneto/ 
ionopause obstacle shapes typical of terrestrial planets. The 
overall objective was the erablemer.t of rational modeling 
studies to be performed on the global so.! ar wind interaction 
problem for the Earth and other terrestrial planets. This would 
involve the development and utilization of a multi-f acited 
ccmputaticnal procedure that embodies the theoretical model 


and incorporates advanced numerical methods so as to enable more 
general and detailed studies to be performed than heretofore 
possible. Comparisons of predictions from the model were antic- 
ipated both with existing theories where possible and with 
observational results where available. Additionally, the model 
was to be employed in several collaborative efforts with other 
space scientists in which predictions from the model would both 
augment other theoretical analysis and assist in interpreting 
results and comparisons with observations. 

3. SUMMARY OF IMPORTANT RESEARCH RESULTS 

The important results obtained under this contract are 
as follows. 

3.1 Development and Verification of the Basic 
Procedure for Predicting ths- Global Plasma and 
Field Properties of the Three-Dimensional 
Solar Wind Flow Past Terrestrial Planet 
Magneto/Ionopause Obstacle Shapes 

The basic theoretical model for predicting the global 
interaction of the solar wind with terrestrial planet magneto/ 
ionospheres involves solution of the continuum partial differ- 
ential equations of magnetohydrodynamics of a perfect gas 
having infinite electrical conductivity and zero viscosity and 
thermal conductivity which express the conservation of mass, 
momentum, energy, and magnetic field. The computational 
method developed to solve these equations is described in 
detail in ref. 1. Briefly, two separate finite-difference 
procedures are employed; one to treat the flow domain encom- 
passing the bow shock and obstacle back to the terminator, 
and another to treat the flow domain downstream of the termin- 
ator. The reasons for doing this are discussed in ref. 1. 

The result is the achievement of a computational efficiency 
not possible by any other existing means. 
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One of ^--he fundamental tasks of this part of the study 
was verification of the accuracy und range of validity of the 
computational procedure. This verification was accomplished 
through comparisons with other theoretical methods, and also 
by exercising the procedure on computations for obstacle shapes 
and solar wind flow conditions that span the entire range of 
interest of applications to terrestrial planets within the 
solar system. These results are reported in ref. 1 and 
demonstrate both the accuracy and robustness of the procedures . 

Finally, one of the most significant achievements of the 
program v;as the incorporation of the entire computational 
procedure embodying the theoretical method into a modular, 
completely automated, user-friendly code. This permits any 
general user access to the results of the model at very reason- 
able computational cost. This feature of the effort is demon- 
strated in ref. 1 where examples are providr3d of detailed maps 
of the global nteraction region for all important plasma and 
field variables. 

3.2 Further Development and Evaluation of the 
Predictive Model by Detailed Comparisons of 
Plasma and Field Properties with 
Observational Data 

The advanced computational model that was developed 
above was further refined and generalized. These developments 
related to both obstacle shape geometry and oncoming solar 
wind conditions and are detailed below. 

A new family of obstacle shapes for solar wind/iono- 
spheric interactions were determined which take account of 
gravitational variation in the scale height. These results 
are reported in ref. 2. A catalog of solar wind flows past 
these new shapes were obtained and results are presented in 
ref. 8. This catalog proved instrumental to a wide number of 
investigators engaged in related solar wind studies. 
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A significant refinement included in the computational 
model under this task has been the allowance of an arbitrary 
oncoming direction of the interplanetary solar wind. This 
refinement is detailed in ref. 2. This feature has been 
instrumental in showing that the asymmetry of the bow shock 
at the Earth previously observed and attributed to magnetic 
effects is either non-existent or at a level below that 
measurable by current spacecraft instrumentation. These and 
extensive related results employing the computational proce- 
cures which further demonstrate the accuracy of the predic- 
tive model are reported in refs. 4 and 5. 

Finally, a provision for the determination of time 
histories of plasma and field properties along an arbitrary 
spacecraft trajectory has been incorporated into the computa- 
tional model. This feature has proved invaluable providing 
theoretical results for direct comparisons with spacecraft 
observations. Details of this provision and results of time 
history comparisons with plasma and field data- from the Pioneer 
Venus spacecraft are provided in ref. 2. In particular, compar 
isons with magnetic field observations for several orbits at 
quiet-time conditions are given and display very good agreement 
Because determination of the frozen magnetic field is the last 
step in the computational model, its accurate prediction serves 
as a critical test of the correctness and accuracy of the 
theoretical and computational models. 

3.3 Collaborative Efforts Employing 
Present Predictive Model 

A number of collaborative efforts have been undertaken 
with other space scientists in which predictions from the 
present model have been employed to interprete observational 
results as well as to augment other theoretical analyses. The 
most extensive of these efforts has been with Dr. James A. 
Slavin and Prof. R, E. Holzer at the Institute of Geophysics 


and Planetary Physics, UCLA. Results from the current pre- 
dictive model provided the theoretical basis of an examination 
of solar wind flows past all the terrestrial planets. This 
was applied, first, to modeling the mean bow shock s^'ape and 
its position and is reported in ref. 4; and then to interpreting 
these comparisons in light of both observations and previous 
analyses (ref. 5). The results from this study reported in 
refs. 4 and 5 represent the most detailed and comprehensive 
study to date of the global interaction phenomena of the solar 
wind with terrestrial planets, from the aspect of observations 
of the bow shock position and shape. They both confirm and 
verify the model's applicability to all terrestrial planets. 

Another collaborative effort employing the current pre- 
dictiv.^ model was carried out with John D. Mihalov of the 
Space Sciences Division at Ames Research Center. In that 
effort, which formed the basis of Mr. Mihalov's dissertation 
for Engineer's Degree from Stanford University, the predictive 
model was employed in a series of calculations. involving time 
histories of detailed plasma and magnetic field properties in 
the ionosheath region of Venus, and comparisons thereof with 
observations from the Pioneer-Venus orbiter. These results 
are reported in ref. 6 and also indicate the correctness of the 
present model. 

Finally, additional efforts have been underway with 
Ms. Nancy Crooker of UCLA regarding ISEE observations. Dr. 
William Knudsen of Lockheed Palo Alto Research Labs regarding 
Pioneer-Venus velocity potential analyzer observations and 
interpretations of the Venusian ionopause, and Dr. Janet 
Lu.hmann of UCLA regarding magnetic field draping observations 
at Venus. 
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Abstract . A coaputational nodal has been 
developed for the detemlnation of the gaadynaaic 
and aagnetic field propertlee of the aolar wind 
flow around a nagnetlc planet, such as the earth, 
or a nonaagnetlc planet, auch aa Venua. The 
procedurea are baaed on an eatablished aingle- 
fluid, steady, dlasipationless, nagnetogasdynanic 
model and arc appropriate for the calculation of 
axlaynnecric, supersonic, auper-Alfvdnlc aolar 
wind flow past a planetary r.'gneto/lonosphere. 
Sanple results are reported ior a variety of 
aolar wind and planetary conditions. Sonc of 
these are new applications; others are Included 
to show that the new procedures produce the aane 
results as previous procedures when applied to 
the same conditions. The new aethods are 
completely automated ard much more efficient and 
versatile than those employed heretofore. 

Introduction 

The magnetogasdynaaic model for the inter- 
action of the solar wind and a planetary magneto- 
sphere, or ionosphere if the planet is nonmag- 
netic, has been found to be of great utility in 
the prediction and interpretation of observations 
in space. The development of this model has a 
long history steaming from the work of Sprelter 
and Jones [1963], Dryer and Faye-Feteraen [1964, 
1966], Sprelter and SuMers [1966], Alkane [1967], 
Dryer and Heckman [1967], Sprelter et al. [1968, 
1970], Alkane and Webster [1970], Sprelter and 
Alksoe [1969, 1970], and Sprelter and Rizzl 
[1974], Although the usefulness of the original 
results is well established, there is auch to be 
desired. Results have been calculated for only 
a limited set of solsr wind snd planetary condi- 
tions. Moreover, the original solutions bordered 
on what was barely possible at tb':; time, required 
considerable hand computation, and were extremely 
laborious. Improvements in both numerical 
methods and computer capabilities have now 
rendered such procedures thoroughly obsolete, 
both from the standpoint of efficiency of calcu- 
lation and generality of application. 

To remedy this situation, a new, efficient 
procedure has been developed that provides the 
nui^ leal solution for the straamllnes, magnetic 
field lines, and contours of density, speed, 
temperature, and magnetic field intensity 
similar to those faailliar from the original 
literature. The basic theoretical model employed 
is that described by Spre. ter, Alksne, and their 
colleagues in the work cited above. The present 
solution mskes full use of recent advances in 
nonlinear computational fluid dynamics. These 
sdvancea Include an algorithm of Beam and Warming 

Copyright 1980 by the American Cmopbysical Union. 


[1976] to determine the gf' .t. lynamlc flow near the 
nose of rhe magneto/ionopause. The ramainder of 
the flow field is calculated by ualsg a method 
dev^zloped by Kutler ct al. [1973] and by Chaussce 
et al. [1975]. The magnetic field is also calcu- 
lated quite dlffetcntly than it was previously, 
being ^aed on the decomposition of Alkane and 
Webstar [1970]. This paper presents an outline 
of the principal features of the analysis and a 
sample of the results. Complete details of the 
computational procedures are provided by Stahara 
et al. [1977]. 


Mathematical Model 


Tl,e fundamental assumption of the present 
work and all of the geophysical work cited above 
is that the average bulk properties of solar 
wind flow around a planetary magneto/ionosphere 
can be described adequately by solutions of the 
continuum equations of magnctogancynamics of a 
perfect gas having infinite electrical conduc- 
tivity and zero viscosity and thermal conduc- 
tivity. These equations, which express the 
conservation of mass, momentum, energy, and 
aagnetic field, are given by the following 
expressions : 


dt 
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where 
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and the et,uation of state of a perfect gas is 
given by 
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In accordance vlch atandard uaage the ayabola b, 
p, V, T, e ■ C^T, and h • e^/, ■ C refer tc 
the denalty, preaaure, velocity, te( stature, 
internal aaergy, and enthalpy; and C refer to 
the apecific hcata at cr >atant voluac and prea- 
aure; R ■ (C “C^)u * 8.31»i0^ *rta/| *K la the 
univeraal conatan'. , and u repreaenta the 
■van aolecular weight nondlaenalonallzed ao that 
. ■ 16 for atoalc oxygen. For fully ionized 
hydrogen, ;> la thua 1/'’. The aagnetic field | 
and the Poyntlng vector S for the flux of electro- 
^gnetlc energy are in Cauaaian unite. The 
gravitational potential • and acceleration | are 
aasuaed to be due to ■laaive fixed bodlea ao that 
their tlae derlvatlvea are zero. Becauae of the 
oalaalon of dlaalpatlve tena in theae equationa, 
aurfacea of dlacontinulty aay develop in the 
aolutlon, acroaa which the fluid and aagnetic 
proper tiea change abruptly, but in auch a way 
chat aaaa, aioaentua, magnetic flux, and energy 
are conaerved. Theae are approximatlona to 
comparatively chin layera acroaa which alallar 
changes in the fluid and magnetic properties 
occur in the corresponding theory of a dissi- 
pative gas and correspond ;'hyaically to the bow 
wave, magncto/lonosphere boundary, and possible 
other thin regions of rapidly changing pro- sreies. 

The shape ot the magneto/ ionospheric obstacle, 
assumed to be axlsy»etric about a line through 
the planet center extending parallel to the 
free-streaa solar wind velocity vector, may oe 
specified in either of two ways: by a numerical 
table of coordinates or by use of simplified 
models from the cited references, which enable 
the calculation of approximate bagneto/lonopause 
shapes upon specification of a small number of 
solar wind and planetary magnetic field or iono- 
spheric parameters. 

For he earth or other magnetic planets having 
.1 dipole field this shape is obtained by rotating 
the equatorial trace of the three-dimensional 


magnetosphere surface determined by balancing 
the Newtonian pressure p ■ Kc cos^^ against 
the magnetic pressure due to !f*°tlmes She tangen- 
tial component of the dipole field. Here, K and 
f are constants usually taken to be unity (see 
Spreiter [1976] for a discussion), c and V are 
the density and velocity of the solar wind, the 
subscript * refers to conditions in the incident 
stream upstream of the bow wave, and is the 
angle between y and the normal to the magneto- 
pause. These considerations lead to the 
following differential equation for the geocen- 
tric distance to the magnetopause r as a func- 
tion of geougnetic longitude t, measured in 
such a way that « is equal to ir/2 at the magneto- 
sphere nose and increases to 3r/2 as one moves 
around the flank of i.he aagnetosphere to the 
remote tail [Bearn, I960]: 


df 
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r/2 < ♦ *. 3n/2 

In this equatlou, r^ • ■ 

D - a (f^B^ /2wltp is the geocentric 
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distance to the magnetoapherc noae, where a 
refera to the radius of the earth, and B^q Is the 
average intenaity of the field at •the geomagnetic 
equator. 

For a nonmagnetic planet having a eufficlcntly 
dense ionosphere to withhold the solar wind, such 
as Venus appears to be most of th* time, the 
ionopause shape is determined by balancing the 
Newtonian pressure and the ionospheric pressure 
approximated by 
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Pr 
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oeriveo trom by assuming steady condi- 

tions, negligible velocity and magnetic field, 
and a constant lonosoherlc scale height H*kT/mg, 
where k ■ 1.38 x 10~^^ ergs/*K it Boltzmann's 
constant, m is the mean molecular mass, and Pg 
is the pressure at some reference radius r^. 

These considerations lead to the following differ- 
ential equation for the coordinates r^ of the 
Ionopause as a function of angle E from the 
ionopause nose as measured at the pla.ietary 
center [Spreiter et al , 1970): 


cir 
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df 


sin 2e - 2 
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0 1 = 1 ■ 
( 9 ) 


in which • exp [-(r^ - Ro)/H], where Rq is the 
value of at the ionopause nose where £ ■ 0. 

A family of ionopause shapes can be obtained by 
integrating (9) for different values of H/R^- 
Appropriate values of H/R^ for application to 
Venus and Mars appear to be in the range from 
0.01 to 0.30. 

Two important parameters which characterize 
the solar wind flow in the magne toga adynamic 
model are the free-stream values for the Mach 
number M - V/a and the Alfvdn Mach number 

« V/A, In which a • (y is the speedy 

of sound, A is the Alfvdn speed A ■ (B^/4it; , 
and > is rhe ratio of specific heats, usually 
taken to be 5/3 in solar wind applications. 
Because is typically greater than 5 in the 

solar wind, an Important simplification may be 
introduced on the basis that for large the 
magnetic terms in the momentum and energy equa- 
tions (2) and (3) are generally sufficiently 
small compared with the gasdynamic terms that 
they may be disregarded therein. This decouples 
the solution for the fl:iid motion from the 
magnetic field, which may be calculated subse- 
quently by using the induction equation (4) and 
values for y already calculated. 

Although few detailed comparisons have been 
made, other than those of Spreiter and Rizzi 
[1974], which made use of their complete magneto- 
gasdynamic numerical solution for the special 
case of aligned flow, it ^s gene^'ally regarded 
chat Che magnetic field calculetcd in this way 
reproduces the broad characteristics of observa- 
tions. However, the process indicates magnetic 
field strengths that may be Coo large iMcdiately 
outside the magnetopause. Zwan and Wolf [1976] 
have made an Intaraatlng analyala of conditions 
in this legion and concluded chat the acrong 
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aagocLlc flalda would drlvt tba ^laaaa out along 
rhc flald linat and produce a thin ragion of 
depleted plaasa denelty. Crooker and Siacoe 
[1977] have ahown further Chat an obaerved 
anlaotropy of the aagnetoaheath preaaurc occurs 
naturally froa a alrror Inatablllty of the 
resulting flow al'‘:.g the aagnctic flux tubes. 
Neither of these properties nor any representa- 
iirr. cf the long discussed possibility of 
aagnetlc aerglng at Che aagnctopcusc are captured 
pv the present aodcl, but the ready ability to 
geners't our solutions for any specified condi- 
tions should aid in further study of such pheno- 
aena. 


Calculation of the Casoy.'iar.ic Flow Properties 

The gasdynaalc equations for supersonic flow 
past a blunt-r.osed obstacle representative of a 
planetsiy sagneto/ionoaphere are solved aost 
effectively by using one aethod for the nose 
region, where the flow accelerates from subsonic 
to supersonic spieds, and another for Che 
reoaining region, where the flow is purely 
supersonic. In the present work the solution 
for the nose region is calculated bv using a 
new axisymtetr Ic lapliclt urwtesdy Eulei *quation 
solver (L*ff), which def.e.Tsi.nes the steady state 
sclution by a liae-aarc hlng procedure. The 
remainder of the sol'. ■ is calculated by 
using a shock-captu > .'ching procedure (SCI), 
which spatially p *■ nces thi* solution downstream 
as far as require'! Ly solving the steady Euler 
equa .on.s. 

laplicit Unsteady Euler Equation Solution 
for the Nose Region 

The partial differential equations esployed in 
the laplicit code are the unsteady gasdynamic 
Euler equations cerlved from (l)-(3) by deleting 
the magnetf and gravitational teras and special- 
izing for axlsymetric flows. Upon introduction 
cf the generalized independent variable trans- 
formation T • T, C ■ C(T,X,R), theae equations 
may be written as 

a- J). * 


♦ Kh^V + r,^E * C - 0 (10) 

where 
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and J » - Cj,hy is the transformation 

Jacobian, in these equations, T denotes tlae, 
X the axial downstream coordinate, and K Che 


cylindrical radial distance; u and v are the 
velocity components in the X and t directions; 

•t " p/|p(T” 1)J ♦ (u* ♦ v*)/2 and aubacripts 
denote partial derivatives with respect to the 
indicated variable. 

The analyala commences by introducing a compu- 
tational mesh in polar (r.B) coord Inatas such 
that one family of coordinates consists of rays 
from the planetary center spaced at equal incre- 
ments of 6 measured from Che obstacle nose and 
Che ocher of curved lines icteraecting each ray 
so as to divide the portion of it between the 
magneto/ lonopause and the shock wave into a 
fixed niaber of equal aegments. The coordi- 
nate transformation equation is than used to 
as.-: the portion of the X,t,T physical space 
bounded by (1) the bow wave, (2) the downstream 
outflow boundary at 6 > v/2, (3) the obstacle 
surface, and (4) the stajnaclon streamline at 
B • 0 into a rectangle in the C, n, t c:«puta- 
tional space as illustrated in Figure 1. Cene" 
ally, the transformation metrics at «ach time 
step are not known beforehand and must be de'er- 
mlned numerically as part of the solution. 
Integration step size is established by using 
the eigenvalues of thc^Jacobian matrices A and B, 
where A - 3E/3U, B ■ if/3V, and also where 
r - U/J, fe - (CtU -a ExE + CrF/J, and 
F • (oyU ■' nxE + ngF)/J. 

Bou^ary conditions for a properly poacd 
mathematical problem arc that the flow satisfy 
the axisymmetric Xanklnc-Hugonlot ahock relations 
derivable from equation (10) along Itam 1 above, 
be entirely supersonic along itam 2 above, be 
parallel to Items 3 and 4 above, and be symaetric 
about item 4 above. Inlt'.al co^iclons are 
determined by use of an ipproxlmat Ing formula 
for the coordinates of the bow shock wsvc, 
dependent on y, H_, and the shape of the obstacle, 
and by prescribing a Newtonian pressure distri- 
butic 1 on the obstacle. The rMslndcr of the 
initial properties on the obstacle surface can 
then be determined from the conditions that t^e 
flow must be tangential to the surface and 
possess constant entropy. A linear variation 
for the flow properties between tbe bow shock 
and the obstacle is then prescribed. This 
provides the initial flow field, ’/hich is ther. 
integrated in a time asymptotic fashion until 
the steady state solution is obtained. 

The basic numerical algorithm uaed in the IMP 
code was developed by Beam and Harming [1976] 
and is second-order accurate, noniterative and 
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Phymcoi plone Compulotion plane 

Fig. 1. Transformation from physical domain to 
rectangular domain. 
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(b) (e) 

Fig. 2. Co»piiri*on of flow proptrtle* predicted 
by pfettnt Implicit method with other technique* 
and experiment for *iver*onic flow past a sphere; 
M « 4.926, y « 1.4. (a) Stagnation line density, 
(b) Surface preaaure. (c) Shock atandoff 
distance. 


S la a fourth-order amoothlng tat* introduced 
to allmlnate noallnear atabllltlaa that may arise, 
since the use of central differences In the 
spatial directions results In s neutrally stable 
algorithm. 

At the boundaries, modification of the differ- 
encing algorithm to account for the particular 
conditions described above Is accomplished a* 
follows. The obstacle surface flow tangency 
condition Is Incorporated through the use of 
Kentser's [1970] scheme, while at the symmetry 
plane the variables are reflected according 
whether they are odd or even. At the outflow 
boundary, where the flow Is entirely supersonic, 
the dependent variables are determined by extra- 
polation from the adjacent Interior points. For 
the upstream bowndary, formed by the bow shock 
wave, the sharp discontinuity approach of Thomas 
et al. [1972] Is used. The Interior flow field 
bounded by there various boundaries Is treated 
In shock-capturing fashion and therefore allows 
for the correct formation of secondary internal 
shocks should any occur. 

Shock-Cspturini Marching Solution 
for the Downstream Solution 


The shock-capturing technique of Kutlcr ct al. 
[1973] and Chaussee et al. [1975] employed herein 
is based on the set of four equations for steady 
axlsymsetrlc flow obtained from (10) and (11) by 


spatially factored. In particular, the 'delta 
form' with Euler time differencing is employed. 
When applied to (10), the algorlttis assumes the 
form 

(I ♦ ird,A“)(I ♦ iTi b")(c"^* - C") 

- - iT(d « f" * G) (12) 

C n 

where A and B are the Jacobian matrices, 1 Is the 
Identity matrix, £. and 6 are, second-order, 
central dlf ference^operatSrs, - C(ntr), and 

iT Is the Integration step size. 

Equation (12) is solved at the Interior points 
only. It requires two 4x4 block tridiagonal 
inversions at each time step of the Integration. 
The solution proceeds as follows: 

1. Define dU - - u”. 

2. Form the right-hand side of (12), and store 

• B+.1 

result* In the L array. 

3. Apply smoothing - (t/8)S/J. 

4. Define U ■ (I ♦ iTi^B'')tC’, and solve the 

matrix equation (1 + i-td.A"*")U • for C, 

^ ^rv4^1 

storing the result In the array. 

5. Folvc the matrix equation 

(I + 1t6 b“)aC' - for dC. 

6. Obtain the values of U from the 

relation • dC 4- C". 

7. Transfer the contents of U to U**, and 

repeat all steps until satisfactory convergence 
Is attained. 



5" 



Flg. 3. Comparison of implicit and inverse 
methods for shock shape and sonic line location 
ar-d density distribution along bow shock and 
magnetosphere boundary for > 8, y > 5/3 flow 
past the rotated equatorial trace of the magneto- 
pause. (a) Shock shape and sonic line location, 
(b) Deaslty dlstributloa. 
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•attlBt ^ ••rlvatlvaa to s«ro. Tb« fourth 
of thla Mt of •^uationa, which roprooonta 
eooaarpatloB of tot*l aaargy pa , can ba 
latagratad to obtain tha ralatlra ^, * *( '*' p/o " 
conat for tha total anthalpy par unit naaa. 

Tha cowputatlonal aaah la daflnad by llnaa of 
conatani X and (r-r. )/(» -E. ) , whara and 
arc functions of X that oaacrlba tha radial 
cylindrical coordlnatas of tha naRnato/lonopauac 
a^ bow ahock wava at tha aana X aa tha flald 
point (X.R) lllustratad In Flgura 1. Tha thrae 
rcBalnlng partla. dlffarantlal aquations for 
conservation of naac and of axial and radial 
aoRcntun arc than tranafomad to a ractangular 
ccaputatlonal apace by C ■ X, n ■ <R-R^)/ (R-R^) 
to obtain 



d «« 


Fig. 5. Magnatopsuae praaaura coafflclanta for 
tha principal naridian nagnatopauaa ahapas ahowr. 
in Flgura 4. 


and h ■ y/fy-l) (p/p) to datanlna tha following 
quadratic aquation for u: 


If- f ♦ G - 0 (13) 

where 
E - E 

c-c— V 

R^ - 

( 14 ) 

The finite difference counterpart of (13) Is 
Integrated with raapcct to the hyperbolic coordi- 
nate ( to yield values for the conservative 
variable E. Subsequent to each Integration step 
the physical flow variables p, p, u, and v wust 
be decoded fron the cowponents m E- This 
neccssltstea the solution of four slaultancous, 
nonlinear equations and la facilitated by using 
the relations v ■ e.,/e^. p - Cj - e^u, p ■ Cj/u, 

2 2 

liaen'tfl snooe 

Free turiKt 

I 6~ 

Efnbeooe: 
pocke' imcMea 

1 6 lurface ani») 

14. 

12 
“ lO- 

08 


0 


I/O 

Fig. 4. Bow atw>ck and Mbaddad ahock locations 
for aolar wind flow with ■ 5, y ■ 5/3 paat the 
rotated principal ■arldlaa'of the aagnatopauae . 




Two roots exist; one corraaponda to aubsonlc flov 
and la diacardad, since u la always supersonic 
In the present application, while the other cerre- 



Flg. 6. Variation of shock standoff distance 

with oncowlng Hach mabar and ratio of specific 

heats for various aagnato/ionopausa traces at 

dataTBlnad by tha present Implicit procedure. 

(a) Magnetopaut • equatorial trace, (b) lonopause 

trace Hl/Ro " 0.01. (c) lonopause trace 

-■/R - 0.5. 
o 
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1/0 


Fig. 7. Shock ahopoo for various auparaonlc 
flova past tba rotatad aquatorlal traca of the 
■agnacopauaa; co«biD*d near (blunt body) and 
far (urchlng) aolutlsna. 

aponds to auparaonlc flow and glvaa the daalrad 
aolutlon. 

Since only the bow ahock wave la traated aa a 
sharp discontinuity and any others that aay be 
present are 'captured' by the difference algo- 
rlcic selection of the approprla.e finite 
dlfferjuce schaae to advance the calculation In 
the i direction la of priae laportance. Aa In 
the blunt body aerodynaalc analyala of Kutler et 
al. [1973] and Chauaaee et al. [1975], the nvaer- 
Ical Integration of (13) la accoapllaed by using 
the finite difference predictor-corrector schrae 
of HacCoraack [1969] the aost efficient aecond- 
order algor Itta for anock-capturlng calculatlona. 

Calculation of the Streaallnea 

The streaallnea are deteralned by Integrating 
trajectorlea through the known velocity field, as 
this procedure was found to be aore accurate than 
the alternative aasa flow calculation. The calcu- 
lation of a particular straaallne la Initiated 
at the bow shock, where Its slope la calculated 
with an exact gaadynaalc relation contained lapll- 
cltly In both the blunt body (Off) and aarchlng 
(SCT) aolutlon, and continued by atepvlse Inte- 
gration In X using a aodlfled third-order Euler 
predictor-corrector aethod. Bivariate linear 
Interpolation froa the flow field grid points la 
eaployed to obtain the velocity coaponents (u,v) 
required at the stepwise points along the streaa- 
llne trajectory. 

Calculation of the Magnetic Field 

With the flow properties Imown froa the gaady- 
naalc solution, calculation of tl>e aagnetlc field 
B proceeds by Integrctlng the steady state 
counterpart of (4). coMonly Intarpratec' aa indl- 



• '• 


Fig. 8. Shock shape for ■ 8, y - 5/3 flow 
past an lonopause shape with B/B >0.1: coabined 
near (blunt body) and far (aarchug) field eolu- 
tlons. 



flow past tba rotatad equatorial traca of the 
aagnatopausa : coabined blunt body and aarchlng 
flow field. 

eating that the field lines aova with the fluid. 
This lands to a straightforward calculation In 
which the vector dlatance froa each point on an 
arbitrarily selected field line to Ita corre- 
sponding point on an adjacent field line In the 
downatreaa direction la datemlned by nuaerlcally 
Integrating fV dt over a fixed tlae Interval 6t. 
Once the coordinates of the field lines are 
deteralned, the aagnetlc field at any point aay 
be calculated froa the relation 

“ pdf /(P-|df_|)i where At la the vector 
length of a aaall aiaaent of a flux tube. 

Such a procedure la valid generally, but Its 
use In tt.e present calculations la confined to 
only one coaponent of the aagnetlc field. The 
other coaponents are deteralned by use of the 
Alkane and Webster [1970] docoaposltlon dependent 
upon the axlsyaaetric properties of the gaadynaalc 
solution and the linearity of the aagnetlc equa- 
tions. 

Calculation of the Contour lines 

Contours are calculated for condlaensionallted 
velocity IY|/V^. density p/p_, and aagnetlc field 


I.. 

■ ’I.C 



Fig. 9b. Velocity asp for M_ > 8.0, y > 5/3 flow 
paat the rotatad equatorial trace of the aagneto- 
pauee: coabined blunt body and aarchlng flow 
field. 
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• oBponenti (•/*.) by application of 

t trarchinc procedure. After a value for the 
;«itour line if epecifiad. the shock boundary is 
•firrhad (or intervals vhich bracket the selected 
T«..i After one such point is located by inter* 
the reasinder of the contour is deter* 
'wslking' around the contour, searching 
. « * step for tne interval and then interpo* 

.’.i tc find the point throush which the contour 
, pasaea. This is repeated until a boundary 
: IS reached. The closed contours are 
in a similar manner. Linear interpolation 
.1 ;.sed throughout the process. The coordinates 
- :■)( contour lines are output as listings, pen 
•.;s. or both. 

Results and Discussion 

Tc verify the correctnets of the new procedures 
s'.c to demonstrate tneir capability for calcula* 
luig solar wind flows for a variety of conditions, 
a large number of test cases have been run, and 
t sacplc 'iiave beer, evaluated by comparison with 
;rrv;ouslv available theoretical and experimental 
rts^lts Figure 2 presents such a comparison for 
s.;trsor.iw flow of air (y>1.4) past a sphere at 
• -.92b. The variation of density along the 
stagnation stream'' ine is provided in Figure 2a, 
vi.-.'.i the variations of surface pressure p. and 
sv.s Standoff distance L with angular distance 
true the nose are given in Figures 2b and 2c. In 
:ht>i and similar comparisons the results provided 
:. thr new procedures are in essentially perfect 
agrrement with those of previous numerical solu- 
tions anj experiment. 

Figure 3 exhibits a comparison of results 
predicted by tne present method with those 
criginally calculated by Spreitcr et al. 11966, 

•vbs for the same solar wind and magnetosphere 
.'onditions, using the inverse method of Inouye 
an; Lomax [1962] for the blunt nose region and 
the method of characteristics for the downstream 
iupersonic region. Figure 3a displays the loca- 
tions of the bow wave and sonic line for flow- 
past the equatorial trace of the magnetopause 
for • 8 and > “ 5/3. The density distribution 
along the magnetosphere boundary and along the 


s. 



1 . t. -i; ; -a. -I. li -s. 


Fig. 9c. Temperature map for M ” 8.0, y ■ 5/3 
flow past the rotated equatorial trace of the 
magnetopause: combined blunt body and marching 
flow field. 
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” 1. t -i. , .•<. -t. -a. -j. 

Fig. 9d. Density map for » 8.0, y ■ 5/3 flow 
past the rotated equatorial trace of the magneto- 
pause. combined blunt body and marching flow- 
field. 

shock wsve are given in Flgu.'* 3b. Essentially 
perfect agreement is obtained between the present 
implicit technique and the inverse method. 
Similarly perfect agreement has been obtained 
with all of the results presented previously in 
the cited references by Sprelter and his 
colleagues. The present method is, however, both 
more efficient and more versatile. 

Figures 9 and 5 prese*-* ne-. results that illus- 
trate the geometric flexibility of the present 
solution and its ability to capture embedded 
shock wsvrs as well as the bow shock, a feature 
which the inverse method cannot duplicate. This 
case is for • 5 and y ■ 5/3 flow around the 
axisymmetric shape generated by rotating the 
principal meridian of an Indented magnetopause 
cbout its axis and may be compared with an experi- 
mental test by Sprelter et al. 11768]. This 
particular profile shape, derived by Spreiter 
and Briggs [1962] from the Beard (I960] approxi- 
mation to the classical Chapman and Ferraro 11931' 
theory, contains s pronounced dent with a concave 
corner in the vicinity of a magnetic neutral 
point. Sprelter and Swmers (1967] argued that 
the presence of an embedded shock wave would make 
it impossible for the magnetosheath plasma to 
follow such an Inde nted contour and that there 
would form instead ,■ free surface approximating, 
in the meridian plans . a tangent line across the 
indentation and capping an embedded cusp region. 
Figure 4 presents results for both the Indented 
and the free surfaces, denoted by solid and dashed 
lines. For the indented surface a shock wave 
is located on the body at approximately 8 • BC* , 
just as in the experiment, while for the free 
surface the pressure coefficient displays an 
approximately constant value as anticipated. 
Finally, we note that the calculation of a 
supersonic flow with an embedded shock and subse- 
quent subsonic pocket provides a severe test for 
any blunt body procedure. The ability of the 
present code to provide convergent results for 
such a flow demonstrates the capability for 
further extension and application to more general- 
ised profiles than was heretofore possible. 

As a final illustration of the range of condi- 
tions for which the implicit code has been tested. 
Figure 6 s«^marlses results for the shock standoff 
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distanca provldad by Cha laadynaalc calculations 
for thraa valuat for y for thrsa dlffarant 
aagaaco/lonopausa shapaa and a ranpe of valuaa 
for froa 3 to 25. Thaaa condltlona span cha 
ranga of Intarast of gaoaatrlcal and solar wind 
conditions for which cha coaputar prograas daval- 
opad haraln would ba nomally applied. 

In ordar to confltw tha ability of tba narch- 
Ing coda to continue the solution to aoae arbi- 
trarily spaclfiad downstraaa location, solutions 
have bean carried far downatraan for a wide 
variety of caaas typical of solar wind Inter- 
actions. Figure 7 exhibits tba location of tha 
bow wave downstraaa froa tha nose region to 
x/D ■ -11 for flows with y • 5/3 and • 3, 8, 
and 25 past tha rotated equatorial trace of cha 
aMgnatopausa. For this calculation, starting 
rondltiona for Cha aarchlng code (SCI) are provi- 
ded by Che blunt body (IMP) coda on tha line 
x/0 • 0.0, which la the usual location at which 
Che two solutions are Joined by the coaputer 
prograa. Tha aarchlng code than dateralnes the 
remainder of cha flow field back to tha speci- 
fied downstrcaa location. Since for the shapes 
considerec herein the downatraas flows arc quite 
ssooth, Che aarchlng calculation is very effi- 
cient (less than 30s, CDC 7600, 0PT«2 compiler). 
Similar results are presented in Figure 8 for 
M *8 and y>5/3 flow past an lonopause shape with 

5 . 




Fig. 10. Contours and field line locations of 
the nagnetic field components (B/B and (B/B^)„ 
in the plane of aagnatlc syMCtry tor “ §.0 
and Y ■ 5/3 flow past the rotated equatorial 
trace of cha magnetopauae . 


n7llg*0.1. Thosa rastilts have bean carried M 
x/ll°"-20 In kaaplag with tha obsarwatlon that 
for Venus and Mara la only slightly greater than 
tba planetary radius, whereas for tha earth, D la 
of tha ordar of 10 earth radii. 

To illuatraf.e tha capability of tha preeant 
procedures to detarmlna streamlines, contour maps 
of various flow properties, and magnetic field 
lines and contours, as well as to dsBonacrata tba 
sutomatsd plotting capability for displaying these 
results. Figures 9 and 10 have bean prepared. 
Figure 9 Illustrates tha coeputar-ganak'stad 
streamline locations and contour naps of velocity 
ratio V/V_, tamparatura ratio T/T^, a,;^ density 
ratio d/p* for the complete near- and far-fiald 
flow about Che equatorial traca of cha magneto- 
sphere for M_ ■ 8 and y>5/3. Based on this gasdy- 
namic solution. Figure 10 exhibits the corre- 
sponding results for the magnetic field (B/B_), 
and (B/B )„ In the plane of nagnetic ajsMstry for 
paralTcl and perpendicular to y . In addition 
CO demonstrating the overall smoothness of the 
computed results these two figures Illustrate Che 
ability of the present techniques to provide the 
completely automated production of report quality 
plots of both gasdynamic and magnetic field 
properties for solar wind flows past axlsyvctrlc 
magneto- ionopause shapes. 

Concluding Remarks 

Advanced numerical techniques have been 
applied to produce an efficient operational 
computer solution for the well-established dissi- 
pationless magnecogasdynamic model for axisyMc- 
tric, supersonic, super-AIfvdnic solar wind flow 
past either magnetic or nonmagnetic terrestrial 
planets. The solution consists of Che following 
assemblage of computer codes; (1) blutic body code, 
to detemine the gasdynamic solution near the 
magneto/ionopauae nose, (2) marching code, to 
determine the gasdynamic solution downstream of 
the magneco/lonopaust nose, (3) streamline code, 
to determine coordinates of flow flsld stream- 
lines, (i) magnetic field code, to uetermine 
frosen-ln magnetic field, (5) contour code, to 
determinu coordinates of contour lines of flow 
and magnetic field properties, and (6) plotting 
code, to plot selected flow and magnetic field 
results. 

Comparisons are reported which demonstrate 
the accuracy of the present techniques by compar- 
ison with previously established theoretical 
methods and with experimental data. New results 
are presented for a variety of solar wind flows 
which Illustrate the flexibility and generality of 
the methods. 
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Abstract . Advanced coaputational procedures 
are applied to an inprovcd model of solar wind 
flow past Venus to calculate the locations of the 
lonopause and bow wave and the properties of the 
<:ovlng lonoshcath plasma In the Intervening 
region. The theoretical method Is based on a 
tingle-fluid, steady, dlsslpatioikless, magneto- 
bvdrodynamlc continuum isodel and Is appropriate 
'cr the calc.ilatlon of axlsyms^etrlc supersonic, 
tuper-Alfvjnic solar wind flow past a nonmagnetic 
plar.et possessing a sufficiently dense Ionosphere 

stand off the flowing plasma above the subsolar 
ocir.t and elsewhere. Determination of tine 
-.iftorles of plasma and magnetic field properties 
i.or.g an arbitrary spacecraft trajectory and 
provision for an arbitrary oncoming direction of 
:r.c Interplanetary solar wind have been Incor- 
tcrated in the model. An outline Is provided of 
:nt underlying theory and computational procedures, 
and sample comparisons of the results art pre- 
sented with observations from the Vloneer Venus 
crbiter. 

Introduction 

The magnetohydrodynanlc model of Spreiter et al. 
!ir7C] for solar wind flow around the Venus Iono- 
sphere, stemming from the closely related earlier 
studies for the earth and Its magnetosphere 
[Spreiter, 1976; Spreiter and Jones, 1963; 

Spreiter and Rizzl, 1974; Spreiter et al., 1966, 
196S, 1970; Dryer and Faye-Peterson, 1966; Dryer 
and Heckman, 1967], provides a theoretical basis 
for the understanding and Interpretation of 
observations from the viewpoint of a fluid rather 
than a particle description of the flow. The 
basic pattern for the associated calculations 
using this model was established by Spreiter et 
al. 11966] for the case of the earth, or any other 
planet with a dominant dipole magnetic field; and 
by Spreiter et al. [1970] for the case of Venus 
and Mars, or any other planet with no significant 
magnetic field but a sufficiently dense Ionosphere 
Co stand off the solar wind above the planetary 
surface. This consists of (1) representing the 
magnetopause or lonopause by a magnetohydrodynamlc 
tangential discontinuity. (2) solving for the 
snape of the magnetospherlc or ionospheric 
otstacle using simplified approximations for the 
pressure of the deflected solar wind flow at the 
obstacle boundary and for the planetary Mgnetlc 
and atmospheric prope-.les, >.3) solving for the 
location of the bow wave and for the density, 
velocity, pressure, and temperature throughout 
the flow as a gasdynamlc problem, disregarding 
the electromagnetic forces because of the charac- 
teristically high values of the free-streas 

Copyright 1960 by the American Geophysical Union. 


Alfv8n Mach number, and flually (4) calculating 
the magnetic field as a frosan-ln convected field 
assuming that the other fluid properties are 
provided sufficiently accurately by the gas- 
dynamlc solution. 

Prior to the work reported recently by Stahara 
et al. [1977], the utility of this model was 
severely restricted because results were avail- 
able only In the form of plots In archival 
Journals for a limited set of solar wind and 
planetary conditions. Results for Intermediate 
points or for other conditions bad tu be deter- 
mined by Interpolation. Stahaia et al. [1977] 
improved on the basic model In many significant 
ways. New and more effective algorithms were 
Introduced to perform the calculations, enabling 
the determination of results for a wide range of 
specific conditions typical of solar wind/ 
terrestrial planet Interactions for both magnetic 
and nonmagnetic planets. 

In the current work suMarlced here, that model 
has been extended and generalised in several 
important directions. These Include (1) extension 
of the gasdynemlc capability for treating free- 
stream Mach numbers as low as M^ > 2, (2) devel- 
opment of a new family of lonopause shapes which 
Includes the effect of gravitational variation In 
scale height of the Ionosphere, and (3) develop- 
ment of the capability for readily determining 
the plasma properties along an arbitrary space- 
craft trajectory, simultaneously accounting for 
an arbitrary direction for the Incident solar 
wind relative to a chosen coordinate system. 

The Hatheaiatical Model - Formulation 
of the Fluid Representation 

The fundamental assumption underlying the 
present work and that reported In all of the 
references cited above la that many properties 
of solar wind flow around a planetary magneto/ 
Ionosphere can be adequately described by the 
contlnum equations of nagnetohydrodynamlcs for a 
single-component perfect gas having infinite 
electrical conductivity and sero viscosity and 
thermal conductivity. 

Governing Equations 

The magnetohydrodynanlc equations for the 
conservation of mass, momentum, energy, and 
magnetic field ere as follows: 

( B fi 

hv^Vk + p«ik - — 
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flux ara coaaarvad. Thaaa ara approx last Iona to 
coaparatlvaly thin aurfacaa acroaa «hlch slailar 
but continuous changas In tha Auld and Mgnatlc 
propart las occur In tha corraspondlng thaory of a 
dlsBlpatlva gas, and thay corraepond physically 
to tha boM vava, lonoiphara boundary, and possibly 
othar thin ragions of 'apldly changing propartlas. 
Tha consarvatlon of nass. aoMntua, anargy, and 
sag-atlc flux laad to tha following ralatlona 
batwaan tha quantitlas on 'ha two sldas of any 
such discontinuity: 



Icv-J - 0 

O) 

where 

[sy • V* ♦ (p ♦ B*^/8t)n - B B /4»] • 0 
n n* c 

(8) 

«1--^- ^ (v'-Vt®ij 

(1=.^ , „ . , . ii] 


and the aquation of state of a perfect gat Is 
given by 

•n 8n^ J 

(9) 

p ■ 16) 

M 

!B • V* - B • V 1 ■ 0 
■ - 1 n n -t 

(10) 


In thase equations and those to follow, the 
s\TEbols c, p, V, T, e • CyT, and h “ CpT refer 
to the density, pressure, velocity, temperslure , 
Internal energy, and enthalpy, and Cy and Cp 
refer to the specific heats at constant volume 
and pressurc_ln the usual way. The symbol 
R • (Cp-Cy)M • 8._3*10’ erg/g*K Is the universal 
gas constant, and M Is the aeon molecular weight 
nondlmenslonalltad so that M ■ 16 for_atoalc 
oxygen. For fully Ionized hydrogen, M Is thus 
1/2. The magnetic field ( and the Poyntlng vector 
$ for the flux of electromagnetic energy are 
expressed In terms of Gaussian units. The 
gravitational potential * and acceleration of 
gravity | are assumeo to be due to massive fixed 
bodies so that their time derivatives are zero. 

C • 6.67*10”®cm Vg s^ Is the gravitational con- 
stant. These equations apply In the region 
exterior to the Ionosphere boundary, as shown In 
Figure 1, and also in a degenerate sense In the 
Ionosphere . 

Because of the omission of dissipative terms 
in these equations, surfaces of discontinuity may 
develop in the solution, across which the fluid 
magnetic proper! 1 's change abruptly, bu * In such 
a way that mass momentum, energy, and magnetic 


Here, (n, t) denote unit vectors normal and 
tangential to the discontinuity surface, where 
q„ represents the local normal velocity of the 
discontinuity surface and ■ v„ - q„ Is the 
fluid normal velocity component relative to the 
normal velocity q,^ of the discontinuity surface. 
The square brackets are used to Indicate the 
difference between the enclosed quantities on 
the two sides of the discontinuity, as In 
[Ql * Q2 ~*)i* ’'here subscripts 1 and 2 refer to 
the conditions on the upstream and downstream 
sides, respectively, of the discontinuity. 
Correspondingly, the average of and Q 2 Is 
represented by 'Q* • (Qj Q2)/2. 

Five classes of discontinuities are described 
by (7)-(10). Of these, only two, the fast shock 
wave and the tangential discontinuity, are of 
significance In the present analysis. The former, 
having properties that satisfy the relations 

V* * C If] >0 Ip] ' 0 [B 1 • 0 
n n 

( 11 ) 

[B )>0 [B-l ■ 0 (.-V*)* .» B'/fi-xl/e') 

• t n n 



Fig. 1. Illustration of overall features of solar wind interaction with nonmagnetic 
terrestrial planets. 
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la appropriate for rapraaanclng the bow wave; and 
the latter, having propartlaa given by 


V* ■ »n “ 0 ’* 0 l?t) * C 

[;) f 0 Ip B^/8*I • 0 


( 12 ) 


is appropriate for repreaenting the ionopauae 
that aeparates the flowing aolar wind and the 
planetary ionoaphere. 

High Alfvdn Mach Nuaber Approximation 


Two important paraaetera that characterlxe 
steady aagnetohydrodynaaic flow paat a atatlonary 
obstacle are the free-streaa Mach number 
M. “ v./a. and the Alfvdn Mach number 
in wnich a • (yp/c)^'^ is the speed of sound and 
A • (B‘/4 »is)^' 2 is the speed of a rotational cr 
Alfvdn wave along the direction of the magnetic 
field, and y ■ fa Che ratio of specific 

heats of the plasma. 

For typical solar wind conditions at Venus, 
both and are high, values of the order of 
10 being representative. As a resclt. Important 
simplifications of the equations may be intro- 
duced. Since tht order of magnitude of the 
momentum flux term in (2) and the kinetic energy 
flux term in (3) is related to that of the 
magnetic terms In the same equations by the 
square of My^, consideration of the magnetic 
terms may be decoupled from the solution of the 
gasdynamlc portions of those equations. The equa- 
tions for the fluid motion thereby reduce approxi- 
mately to those of gasdynamlcs, and the magnetic 
field 8 can be determined subsequently by solving 
the remaining equations using the values for o 
and V so determined. Furthermore, because of the 
Is.'ge value for H^, the terms involving g and t 
can be disregarded in the solar wind flow because 
of the relatively small effects of gravity on the 
flow. This is not true, of course, in the iono- 
sphere, where the comparative stationarlty of the 
gases results in a very nearly hydrostatic 
variation of the pressure with altitude. 

The precise specification of the range of use- 
fulness of the decoupled approximation remains 
elusive in the absence of exact solutions of the 
complete magnetohydrodynamlc equations. The only 
such solutions for the present category of flows 
are those of Sprelter and Rl<:zl [1974] for the 
special case in which B_ is aligned with y^. 

They indicate results virtually identical to those 
of the simpler decoupled theory when is about 
10 or greater and notably differing results as 
is decreased below 5. 

For other magnetic field alignments one can 
only make estimates of the relative magnitudes 
of the divergences of the neglected magnetic 
terms BiBi(/4it -f B^di^i/Sn and compared to the 
terms CVjVj^ + pi^jy and cV|j(v2/2 e + p/e) 
retained in the gasdynamlc analysis and Chen 
speculate about how large the magnetic terms can 
become before they invalidate the solution. In 
regions where the momentum flux and kinetic 
energy terms are large in comparison with the 
pressure and enthalpy terms, the souare of the 
local Alfv^n Mach number M]^ - (ov^/2) > (B^/Bn) 


may ba coosldared to provide a rough measure of 
Che dominance of the gasdynamlc term* in the 
calculation of the fluid motion. .Since this 
quantity is related to Mj[^ by My^ « (o/o^)(v/v^)*. 
(B^/g)*My^, and o/p_, v/v^, and B/B^ are Indepen- 
dent of at large Myy^, it follows Chat My^ 
tends to rasMln large for large My^ except where 
pv* « p_v^ or where B* » B^. These conditions 
are most likely Co be encountered near the 
stagnation point, where v vanishes; but p reaches 
a maximum value comparable to p v^ there, and 
that contravenes the initial premise upon which 
the estimate is predicated. A localised region 
of small My^ does not therefore necessarily signal 
a failure of the approximation. 

As a further consideration in making such 
estimates, it should be noted chat the orienta- 
tion of B^ is very Important in Che determination 
of B^/B^ throughout Che flow. In particular, 
near the stagnation point, b’/B^ » 1 if B. makes 
a substantial angle with v^; whereas B^/B^ 
vanishes ar Che stagnation point when ( is 
parallel to y„, since B is everywhere proportional 
to py [Sprelter at.-f Rizzl, 1974]. For this case 
therefore Mj^ ■ (p./o/M^, end My^ does not vanish 
even at the stagnation point because p/p^ is 
about 4-5 for reprasentatlve conditions. 

A definitive theoretical evaluation of the 
decoupled approximation will have to await the 
development of more advanced methods for the 
solution of the complete magnetohydrodynamlc 
equations. In the meanwhile, conaiderable in- 
sight can be gained by comparing the predictions 
with conditions actually measured in space during 
appropriately steady conditions. The small scale 
of the Venus flow field compared with that of the 
earth's magnetosphere makes the data from Pioneer 
Venus particularly valuable for this purpose; and 
the results presented here show our initial 
efforts to carry out such comparisons. This 
should be considered no more chan a beginning, 
however, as we only have data from two orbits; 
many additional romparlsons should be carried 
out tc obtain a more comprehensive understanding 
of the relation berween the checretlcal and 
observational results. 

Under the conditions outlined above, (2), (3), 
(8), and (9) are thus simplified to the following: 


3^ <"Vk p‘ik> ■ ° 


3 fpv^^ K 3 

3T [— + 


+ e p/p 


] - 0 
J 


[py • V* + p] - 0 
[v* • + pe + p| j « 0 


(13) 

(14) 

(15) 

(16) 


Theje equations, together with (1) and (7), are 
just those of gtsdynamics and mtist be solved for 
the location of the bow wave and the gasdynamlc 
properties y, p, p, and T of the flow field about 
Che ionospheric obstacle. With these properties 
so determined, ,4) and (10) can be uaed to deter- 
mine the magnetic field B. Although all of the 
solutions daacrlbed here are for the ateady state 


7 /IS 


Sprattar and Scahara: Solar Wind Flow Faas Vaoua 



Fig. 2. Illuatratlon ot lonopauaa shapat for atitoa^haraa with varloua (1) conat.wC 
■calc halght* H/I^, and (2) gravitationally varying acala halghta aaaoclatad with a 
constant atnosphcrc and charactarlred by fi/Ro- 


obtained by setting i/3t ■ 0. and v* • v„. l.c., 

• 0, the unsteady equations are prasentsd In 
the foregoing because one of the coapiitac ional 
DCthods used to dctemlne the gaadynanlc solution 
eeploys an unsteady procedure. Integrating In 
time until the steady state solution Is obtained 
asyaptot Ically . 

Det er alnat icn of the Ionosphere Boundar y 

The deteralnatlon of the ionoaphere boundary, 
or lonopause, isakea uaa of the Idealization that 
the lonotpncrlc plasK-i Is sphsrlcaily synoecrlc, 
hydrostatically supporttd. Infinitely electrically 
conducting, devoid of any stagnetlc field, gravi- 
tationally eoun-* to the planet, and Incapable of 
mixing with the aolar wind plasma, as sketchsd 
in Figure 1. This Intsrlot plasms Is censidarad 
CO be teparatad from the flowing aolar wln< 
plasma by a tangential discontinuity across 
which the relations of (12; apply. Because the 
agnetlc pressure la much smallar than the 

Ionized gas pressure p In the lonopause, the 
rtilstlon ip B*/8") • 0 of (12) may be 
simplified to 


Pat. • ^ ®‘'®'>flow 

The variation of the pressure with radial 
distance from the planet center la given accord- 
ing to the hydrostatic aquation by 


(18) 


where p and i art the Ionized gas preaturc and 
density, r la the radial distance measured from 
the center of the planet, and g> ||{ *g,(r,/r)^, 
where aubacrlpt a denotes valuca at the aurface 
of the planet. Combination of (18) and the 
perfect gaa lav given by (6) provides 

p ■ Pr **p [- 1 f ] 


where pj^ ls_tha pressure at some reference radius 
Rg and H “ RT/Mg represent s the local scale 
height of the Ionospheric plasma. 

If H la regarded as constant, l.e.. If 
variations of g and T with r are dlsrcgardeil, 

(19) can ba Integrated directly to obtain 


P “ Fr «*P 



( 20 ) 


Becauac of limited knowledge of the ionospheric 
propertlsa of Venua at the time, the simple 
variation of p with r given by (20) was adopted 
In the Initial study of Spreltar at al. (1970) 
and also In the more recent analysis of Stahara 
et al. [1977] Involving the Initial application 
of advanced computational mtehoda to this problem. 
Data from the Plonear Vsnus spacecraft [Knudsen 
et al., 1979a, b] have provldsd much additional 
knowledge of conditions In the Venus Ionosphere 
and. In particular, have shown .hat the assump- 
tion of an Isothermal (T ■ const) Ionosphere at 
tht heights of concam bars Is quits rsasonable. 

V .a this assumpt on. Including the variation of 


d^i/dr - -eg 


OWOlMAL PAGE B 
OF POOR QUALITY 


Sprtltcr and Scal^ra: Solar Hind Flow Paac Vanua 


77i9 


fra'-lty with height, the following reault for the 
preeaure ia obtained: 


P ■ Pjj «xp 


•s, • - V 


H • r 


(21a) 


V' 


t-rt- 


H • H 

s 


(R^/R,)^ 


(21b) 


iv. oeing the planetary radiua and Hg ■ RT/iigg. 
Equations (20) and (21) provide the two Bodels 
eeployed in this atudy for the ionosphere pres- 
sure variation required in the pressure balance 
relation of (17). 

On the basis of the typically hypersonic 
values for It has been assuoed in all the 
previous analyses cited above that the pressure 
cf the solar wind plasna on the lonopause could 
be represented adequately by the Newtonian 
approximation 

2 : 

p • K. V cos . 

K S 


K 


1 

1 lC->-l)/2] ^ 

'' >-('-l)' 2 xf 


(22) 


where is the angle between the outwaru normal to 
the lonopause and the flow direction of the undis- 
turbed solar wind incident on the bow wave and K 
is a constant usually taken as one but whose 
actual value to provide the correct gasdynamic 
stagnation pressure at the nose of the lonopause 
is as in'* cated. 

The Newtonian approximation stems from gas- 
dynanic analysis of hypersonic flows and is 
appropriate in the range of v from 0' to somewhat 
less than 90°. In view of the close association 
between the magnetic pressure B°/8r and the gas 
pressure p in magne tohydrodynamlcs , it seems 
equally plausible to assume that the corresponding 
relation 


1 T C 

p + B*/8r « K;_^v2ccs“', (23) 

is an appropriate, if not superior, relation for 
the combined magnetic and gas pressure of solar 
wind plasma on a planetary lonopause or magneto- 
pause . 

For the high Mach number flows typical of '>olar 
wind conditions, K approaches 0.881 for 1 * 5/3 
and 0.84<> for > > 2. A discussion of effects of 
nincr constituents, such as ionized helium, in the 
solar wind and of some differences between a fluid 
and collisionless particle representation has been 
given by Sprelter [1976], but the most important 
effect of the introduction of the Newtonian 
approximation is that 'he calculation of the shape 
of the lonopause decnioles from the calculation of 
the properties of the exte-nal flow. In this way 
we arrive at the following equation for the pres- 
sure balance at the lonopause location R^: 

K.^v‘cos‘. • Pj{’(R^) (24) 

where 



depending upon whether or not the gravitation 
variation is included In the scale height. It is 
convenient to choose the reference redlus Rg to be 
the diatence Rq from the center of the planet to 
the nose of the lonopause. Then Pg ■ Pg ■ Kr> v^, 
and " “ 


2 

cos 


i 


A(R^) 


(26) 


at all points along th« lonopause. To proceed, 
cos’ti must be expressed l-i terms of R^ and 6, where 
R^(£) represents the radial coordinate of the 
lonopause and 8 is the angle mesaured at the center 
of the planet with respect to a line that extends 
directly upstream. The appropriate expression is 
[Sprelter et al., 1970] 


2 

cos I, 


i dS J 


(R^d&cosB + dR^sin£) 
dR^ + (Rjde)^ 


(27) 


Substitution of this expression into (26) results 
in the following ordinary differential equation 
for the ordinates of the ionosphere boundary: 




Fig. 3. Comparison of former and present cotqiu- 
tatlonal procedures for determining the gas- 
dynemlc flow properties of solar wind-magneto/ 
lonopause Interactions. 


7720 


Spr«lt*r and Suhara: Solar Wiad Plow Paat Vanua 


ORIGINAL PAGE 18 
OF POOR QUAUTf 


dB 


aln2B - ih - 
2(A - alD^B) 


0 1 0 s. « 


(? 8 ) 


Rcsulta for various ionopausa shai-as obtainad 
by Intagratlrf (28) for dlffarant valuat of H/>c 
using tha cotxatant scala halghc aodal aquation 
(2Sa) wara provldad by Spraltar at al. [1970]. 
Slallar raaults using tha isothaiiaal aodal aqua- 
tion (2Bb) In tha ranga 0.01 i. H/K^ ^ 0.5 ara 
provldad In Figura 2, whara for coMparlson tha 
constant scala halght shapas for corraspondlng 
H/Rq valuas ara also lllustratad. Data froa 
Plo.taar Vanus indicate that tha range rf Interest 
for sppllcatlon r.c Venus appears to be 
0.01 a H/Rg i. 0.10. They also Indicate that tha 
theory provides a good representation of the 
lonopause shape on the daysldc but that conditions 
on the .nightsidc arc both aora coaplax and variable 
than presently accounted for by the theory (see, 
for cxaaple. Brace at al. [1979], Russell ct si. 
[1979b], or Brace et al. (this Issue). To pro *'de 
insight Into soae of the changes to be encountered 
In the theory, sooe results are also Included 
herein ir which the nlghtslde lonopause Is 
arbitrarily flared In and out froa the calculated 
shapes of Figures 1 and 2. 


Calculation of t.ie Casdynaalc 
Flow Properties 

Setcralnatlon of tha gasdynamte flow properties 
Is the most difficult and tlac-consualng portion 
of the numerical solution of tha foregoing equa- 
tions. As sketched In Figure 3, two different 
laathods are used In tha calculations: one for the 

nose region, where both subsonic and supersonic 
flows occur, and another downstreaa of that 
region, where a acre coaputatlonally econoalcal 
procedure can be aaployed, since the flow Is 
everywhere supersonic. In the original analysis 
of Sprelter et al. [1970] an Inverse iteration 
tnethod was used for the nose region, and the 
method of characterlati<*s was employed for the 
remaining supersonic region. Both of these 
procedures are Inferior to more recently developed 
methods and have been superseded In the current 
model. The nose region Is treated u«lng a new 
axlsyxaetrlc Implicit unsteady Euler equation 
solver specifically developed for the present 
application. That procedure determines the 
steady state solution In the nose region by an 
asymptotic tlme-marchlng procedure. The down- 
stream solution Is determined by a shock capturing 
marching proesdure which spatially advances the 
solution downstream as fsr as required by solving 
the stesdy Euler equations. Inasruch as these 
methods are complex and lengthy to describe, only 
an outline of the procedures Is provided here. 


The Sose _Reil on Solution 

The nose region procedure solves the unsteady 
gasdynamlc equations (1), (6), (7), and (13)-(16) 
for axlsymmctrlc flow. The analysis coMcnces by 
mapping the portion of the lonosheath bounded by 
the stagnation streamline, the bow wave, the 
lonopause, and a line through the terminator 
(6 ■ w/2) Into a rectangle. The physical boundary 
conditions are also transferred appropriately. 
Initial flow field conditions are prescribed by 


using a linear varlatlom of the flow properties 
between chose calculated at an aesumsd location 
for the bow wave and Chose calculated at the lono- 
pauae, using e Newtonian pressure distribution 
approximation on the latter boundary. Starting 
with this Initial flow, tha unste.*.dy equations are 
Integrated In a tlme-asymptotlc fashion until Che 
steady state solution Is obtalueu 

The basic numerical algorithm employed was 
developed by Beam and Warming [1976] and Is second- 
order accurate, noniterative, and spatially 
factored. At ttie lonopause and bow shock bound- 
srles a modiricatlon of Che differencing algorithm 
Is Introduced to sccount for the .(articular 
features of the present appllcai Ion. The tangency 
condition at the lonopause is treated u. Ing the 
method of Kentzer [19'70], and he sharp dlscor- 
tinulr.y approach of Thomas er il. [1972] Is used 
at the bow wave. Any secondar ■ shocks which may 
occur In the nose region arc aitomatlcally 
determined via the shock-capturing characteristics 
of the algorithm. 

In our earlier work [Stahara, et al., 1977] It 
was found that for certain lonopcuac shapes which 
have a significant amount of lateral flaring at 
the down-dusk terminator, for exaa^>le, constant 
scale height shapes for H/Rg 2 . 0.5, and/or cases 
involving free-atream Mach numbers a. 3, 
the axial >ior onent of velocity at sosm. points on 
the termina.or plane S * v/2 nay become subsonic. 
Although this has no effect on the nose region 
solver, the downstream solution cannot be obtained 
for these cases, since the march Ing- region solver 
'.-hich determines the solution downstream of this 
starting plane requires supersonic axial veloc- 
ities In order to proceed. In the present effort 
this limitation has been removed by developing the 
capability fur adding an additional portion of the 
flow field, located downstream of the terminator, 
to the blunt body solution as illustrated in 
Flg'..re 4. This effectively generalizes the 
capability of the present procedures to treat a 
wide variety of lonopause shapes. Including all 
of the shapes described by the constant scale 
height and constant temperature models found using 
(25a) .ind (25b), as well as to treat lower free- 
streom Mach numbers down to about M >2. 

The Downstream Region Solution 

The solver for the downctream supersonic region 
is based on the original work of Kutler et al. 
[1972, 1973] and Chaussee et al. [1975] and solves 
for axisyMetrlc flow the steady state gasdynamlc 
equations In conservation law form Indicated by 
(1). (6), (7). and (13)-(16) with the time 
derivatives set to zero. Of these, (14), which 
represents the conservation of energy, can be 
integrated for the steady flow to yield the 
following equation for the total enthalpy: 

h^ " h (u^ * v*}/2 • const (29) 

where h • e '4 p/o • CpT Is the enthalpy per unit 
mass. This relation also holds In the nose region 
In the ll.sit as the steady state flew solution Is 
approached. The computational mesh defined by 
lines of constant X and (R - R^)/(R^ - R^), where 
R, and R^ are functions of X that describe the 
radial cylindrical coordinates of the bow sliock 
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Ionosphere boundary 
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Fig. 4. Illustration of additional flow field segment to the obstacle nose solution 
for determining the gssdynamlc flow properties of solar wlnd-ionopause interactions. 
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wave and ionopause at the same X as the field 
point (X, R), is transformed into a rectangular 
computational grid. The gasdynamic equations are 
similarly transformed and Integrated numerically 
to obtain the solution. Since only the bow shock 
wave is treated as a sharp discontinuity and any 
others that may be present are 'captured' by the 
difference algorithm, selection of the appropriate 
finite difference scheme to advance the calcula- 
tion in the X direction is of prime Importance. 
Following the analysis of Kutler et al. [1972, 
1973] and Chaussee et al. [1975] the method used 
is the finite difference predictor-corrector 
scheme of MacCormack [1969], the mc/st efficient 
second-order algorithm for shock-capturing 
calculations. 

Calculation of the Streamlines 


and the streamline determination is made by step- 
wise integration in X using a modified third-order 
Euler predictor-corrector method. Bivariate 
linear interpolation between the flow field grid 
points is employed to obtain the velocity com- 
ponents (u, v) required at the stepwise points 
along the streamline. 


Calculation of the Magnetic Field 


With the flow properties known from the gas- 
dynamic calculations, the magnetic field B may be 
determined by integrating (4) and (10); or 
alternatively, the following equations may be 
derived from them: 


JL (( 

Dt JJs 


B'dS - 0 



Dt 


- (B-V)y 

P - • “ 


( 32 ) 


Th'f accurate determination of the gasdynamic 
streamlines is crucial, since they provide the 
basis for the frozen magnetic field calculation 
described in the following section. The stream- 
lines are determined by integrating fluid particle 
trajectories through the known velocity field, 
since this procedure was found to be more accurate 
than the alternative mass flow calculation. The 
calculation of a particular streamline is 
Initiated at the point where the streamline 
crosses the bow shock. At that point, exact 
values of the streamline slope dRg/dX are given 
by the following relation Involving M^, Y, and the 
angle between the free-stream direction and the 
tangent to the shock wave: 

dRj (2cotS.)(M^sin^4. - 1) 

— ^ - = = (30) 

dX 2 ■> 

2 + M (Y + 1 - 2sin‘£c) 

which is contained implicitly in both nose and 
marching region gasdynamic solutions. Ac other 
points, the local streamline elope is given by 
the ratio 


dRg/dX - v/u 


( 31 ) 


xn which S is an arbitrary surface moving with the 
fluid and D/Dc is Che material or substantial 
derivative dD/Dt “ 3/3t + (v*V). These equations 
are consonly interpreted as indicating the field 
lines move with Che fluid. For Che steady state 
in which 3/3t • 0 and v„ ■ v„, these equations 
lead to a straightforward calculation in which the 
vector distance from each point on an arbitrarily 
selected field line to its corresponding point on 
an adjacent field line in the downstream direction 
is determined by numerically integrating / y dt 
over a fixed time Interval dt as indicated in 
Figure 5. Once the magnetic field lines are 



Bov ioo:> 


lonopouse 


Fig. 5. Illustration of quantities used for mag- 
netic field line calculation in the plane of 
magnetic sywetry. 
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dataralnad, cha aagoatic flald at any point nay la 
calculatad fro* tha ralatlon 


B 




(33) 


where la the vector length of a aaall aleaent 
of a flux tube. 

Such a procedure la valid generally, but Ita 
uae la confined In the prerenc calculation to only 
the conponent of the aagi.ttic field (B' chat Ilea 
In the plane of augnctlc ayaaetry deflu..d by the 
plane containing Che axla of ayoaetry of the lono- 
pauae and the magnetic field lines upatreaa of the 
bow wave and furthermore la aaaoclated with only 
the component of the interplanetary magnetic field 
''hat Ilea in Che plane of magnetic ayvaetry and la 
perpendicular to y_. The remainder of the magnetic 
t^ld calculation makea uae of a decomposition due 
ro Alkane and Webster [1970] In which Che 
axisynmetrlc properties of Che gasdynanic solution 
and Che linearity of the magnetic field equations 
(4) and (10) or (32) are employed to derive the 
following relationship for the magnetic field Bp 
at any point P: 


?P 





+ e 


B 

B , » 


(34) 


As Illustrated in Figure 6, subscripts ", and n 
refer to contributions associated with the com- 
ponent B," of parallel to y_; the component B«i 
perpendicular to v^ In the plaue that contains the 
point P, the center of the planet, and the vector 
V.; and the component normal to Che latter 
plane, where e„ is a unit vector In r,he latter 
direction. The unit ratios (@p/B,)„ and (Bp/B^)n 
can be calculated directly from the gasdynamic 
solution by the expressions 


\.Kj„ ' c V ! ’ ' 'S> 


(3:.) 



Fig. 6. Illustration of the componr ’ts of the 
three-dimensional magnetic field. 


where Ip la the radial cylindrical coordinate of 
the streaaillne through P, aa Indicated In Figure 6. 

In carrying out tha determination of (Ip/B.)^ 
uslitg (33), values for are determined 

initially at the points ^ere the atraamllnas and 
perpendicular-component flald linea intersect. 

A generalized quadrilateral Interpolation scheme 
fo.'.lowed by a flfth-ordar smoothing la then 
employed to determine cha corresponding values at 
the computational grid points where values for 
d/p. are available for calculation of (Bq/B.)x. 

At the bnw shock, axact fore<jlaa are used for th.' 
magnitude and dlreckion of 


(|.4£!/li4.|)^ - 1 ♦ cot^e(l-hD^) 

- 20 » CSC 9 » cot 6 « coa(9-6) 

(36) 

“ e + sin ^ { [D » cot 9 
» aln(6-6)]/(|if l/lif^l)) 


where 


0^ - 1 - 6{M^sin^9 - 1) 

X (>Mfsin^9 + l)/I(Y+l)^Sin^9) 


cot6 - ^(Y*l)M^/(2(M^sin^9- D) - ij tan9 

a . -^1 


(dXj 


(37) 


and li Is the angle between the free-atrean direc- 
tion end the !jt vector. 

Finally, the resultant magnetic field vector 
can be expressed In terms of components relative 
to any orthogonal (X, Y, Z) coordinate system. 

If the coordinate system Is oriented so chat the 
point P lies In the (X, Y) plane, the magnetic 
components are 


Bv/B “ [('b’/» )„ cos J cos a 
X * * p 

-> (!b'/b ). cos i' sin a )]cos a 
- * « i p n 

B„/B • KIbI/B )„ sin j cos (38) 

y w • * • p 

+ (|b|/B ), sin ii sin a )jcos a 
• * * p n 

B,/B - (B/3 ) sin a 
i “ • n n 

where t Is the local flow angle given by 

1 “ t«n ^ (39) 

and tl.e Intcrplane'ary magnetic !i id angles a 
and a Indicated in Figure 6 art,- defined bv ** 



» • 


r ^ 

-1 

B 

• 

A 

-1 


a - tan 
P 

B 

*M 

■ t«n 

03 
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Fig. 7. llluatration of aun-planct (Sg, y,, Cg) and solar wind (x, y, t) coordinate 
syatcms and the azlnuthal (n) and polar (dp) solar wind anglaa, both shown In a 
positive sense. 



(41) 


Th* generalisation of these results »rt>an the point 
P is at sone arbitrary (X, Y, Z) point not In the 
(X, Y) plane are provided below In the spacecraft 
trajectory section. 


Calculation of the Contour Lines 


In order to provide convenient suasary Inforaa- 
tlon of the detailed plasna and nagnetlc field 
properties datemlned by the above nodel, contour 
and field line waps are deterwlned for the Iwpor- 
tant Held properties. Contours are calculated 
for nondlacnslonallzed velocity |y|/v_, density 
ohm, and aagnetlc field cooponente 

and (B/B_)„. Details of these calcula- 
tions are provided by Stahara at al. [1980]. 

Since the teaperature Is a function of |y|/v^ only 
for a specified H and Y, i.e., 




(42) 



Fig. 8. Illustration of solar wind (x. y, z) and (X, Y, 2) coordinate aystesu and 
the Interplanetary aagnetlc field and nagnetlc field angles (a^t o^) . 


k 
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*■* Mal.O r«ll$7 H/R.»0.10 
STRCOMlINEJS 





(a) 


Mig.c ysi'.b? 



Fig. 9. Calculatcd_lonopauaa and bow shock locations and lonoshaath propartles for 
• 8.0, 1 ■ S/3, H/Rg • 0.10, and altbar parallal or perpandlcular to v_. 

(a) Streanllna nap. (b) Velocity nap,*‘(c) Danslty nap, (d) Tenperatura nap. 

(e) Parallel taagnctlc field nap. (f) Perpandlcular magnetic field nap. 


velocity contours nay also ua considered as tem- 
perature contours with only a relabeling required. 

Solar Ecllptlc/Solar Wind Coordinate 
Trcns format Ion 

To compare rcsul*'s of the present calculations 
with observational data. It la necessary to con- 
sider the appropriate transformations between the 
spacecraft and solar wind coordinate aystens. 

Part of the data required as Input to the theo- 
retical model consists of oncoming Interplanetary 
values of solar wind tenperature, density, and 
velocity and magnetic field vector components. 
Inese are usually reported In a sun-planet or 


solar ecliptic reference frame, whereas the 
natural coordinate system for the theoretical 
model Is one which aligns the axial direction with 
the oncoming solar wind, slnca the gasdynamlc 
calculation Is assumed to be axlsymmetrlc about 
this direction. Thus the Interr' tnetary Input 
data must be transformed to the solar wind system 
to Initiate the theoretical determlnutlon. Once 
the gasdynamlc and magnetic field calculations In 
the solar wind system are complete, those results 
must then bw transformed back to the sun-planet 
system to allow direct comparison with the space- 
craft data obtained as a function of time along 
Its trajectory. 

For the sun-planet system (xg, yg, Sg) the 
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(c) 



Fig. 9. (continued) 


orlgir la placed at the planetary canter, the 
ax. 8 la directed toward the aun, the y, axla 
la directed opposite to the planetary orbital 
action, and the Zg axis la orthogonally northward. 
The direction of the oncoming aolar wind la such 
that the azimuthal angle, including aberration 
effects due to planetary motion, of the solar wind 
velocity vactor In tha planetary orbital plan* la 
•*., and the out-of-orbltal plane polar angle la dp. 
The positive sense of the azimuthal angle la for 
weat-to-eaat flow, and that of the polar angle la 
for north-to-south flow, as Indicated In Figure T. 
Also Illustrated Is the solar wind (x, y, z) 
coordinate system defined by (n, dp). This system 
is somewhat Inconvenient because the direction of 
tslar wind flow Is In the negative x direction; 
therefore the Internal gasdynamlc and magnetic 
field calculations are performed la the (X, Y, Z) 
system Illustrated In Figure 6. 

The coordinates and vactor transformations from 


the ecliptic sun-planat (Xg, yg, Zg) system to the 
(X, Y, Z) solar wind system are given by 


[V 


-cosHcosp 

p 

-sinneos^ 

P 

alnt 

P 


Qxg 

% 

• 

a Inn 

-cosn 

0 


‘Jys 



-cosHsln^^ 

slnnslnt 

P 

cosdp 


r-j 


where (Qj^, Qy, Q^) rapraaants tha Cartesian com- 
ponents of any vector referred to the (X, Y, Z) 
system, and (Qx,i Qjg, Qr_) represents the 
corresponding vecter In the (xg, yg, Zg) system. 
Thus for a transformation of coordinates. 


(Qx- Qy. Qg) - U. Z) 

(Q, . Qy . Q ) - (X,. y,. .,) 


( 44 ) 
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Fig. 9. (conctnucd) 


while for a vector tranaformatlon of, say, the 
magnetic field. 


(Q,. Qy. Q^) - (B,. By, B^) 

• Qy ^ Q, ) ■ (»x • 8y • ) 

s a a ^a a 


(45) 


The inverse transformation from the solar wind to 
the sun-planet systam la given by 


(q ) 

' *al 

-C08^08^ 

P 

slnfl 

-C08H81D4 

P 


Qx' 

Qv i ■ 

' 8 

slnncoi# 

P 

-cosO 

tlnrisln# 

P 


Qt 

N 

O' 

8ln^^ 

P 

0 

COSPp 


*^2 


(46) 


Properties Along a Spacecraft Trajectory 

One of the primary alma of the present effort 
has been to develop the capability to determine 
plasma and magnetic field properties, as predicted 
by the theoretical modal, at locations specified 
along an arbitrary spacecraft trajectory, and In 
such a form as to enabla comparisons to be made 
directly with actual spacecraft data. To 
accompllah this, the following procedure has been 
developed and Implamanted In the current model. 
First, from the known oncoming interplanetary 
condltlona provided In a sun-planet reference 
frame, the azimuthal and polar aolar wind angles 
(n, pp) are employed In (43) to eatabllsh both the 
location (Xp, Yp, Zp) of the trajectory point and 
the Interplanetary magnetic field componenta 
(Bg., By., Bg.) In cha solar wind (X, T, Z) frame. 
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Xonopausc H/K *0.03 


Fig. 10. Illuatration of typical flow field grid dcnalty for gaadynaalc aolution; 
- 3.0, y - 5/3. 


'.n the next etep, in which the axiiysmetr ic gas- 
dvnamic and unit magnetic field calculations are 
carried out, it is convenient to lake advantage of 
the axlsysnetric property of the gaedynamlc flow 
and introduce still another coordinate system 
(X, y', z') related to the (X, Y, Z) system by a 
rotation about the X axis by an angle f such that 
the >■' axis passes through the point P. This 
rotation defines a new coordinate system (x', y', 
z'), where 


fl 


'1 

0 

° 1 

fx] 

0 

cose 

I 

■ in6 1 

Y 

0 

-sine 

cose] 



111 which 


e » tan'^ ^ 

. P 


Thus the (x*. y') plane which contains the X axis 
and the arbitrary point P corresp onds directly tc< 
the plane (X. R)’- (Xp, /Y* + Zj) , in which the 
axisyoBietric gasdynaoic flow properties arc calcu- 
lated. In particular, the valoclty magnitude v, 
density o, and flow angle « at the point P are 
found by bilinear Intarpolatlon through the (X, R) 
flow field grid. The vector velocity in the (X, 

Y, Z) system is then given by the transformation 


- 0 cos8 -alne V sint 


^v^j sln6 cose I 0 J 


and then in the solar ecliptic system by the 
transformation Indicated by 


v -COSliCOSi 

V • slnOcost 

y I p 


V sint 

, *»J I ” 


sinr. -CCS us ins 


.osT, slnHslnS 


v' - A- ♦ Z- 
P P 

z' • 0 


Calculation of the magnetic field at an 
arbitrary point is somewhat more complicated, 
since these components are dependent upon the 
orientation of the incident interplanetary mag- 
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wind dlrectlona; also, various bow shock shapes for different Interplanetary solar 
wind conditions. 




Fig. 12. Cc^parlibu of observed (OPA) and theoretical tiac histories of lonosheath 
plasma properties for Pioneer Venus orbit 6 based on Inbound and outbound Inter- 
planetary solar wind conditions using a gaadynaalc solution for ■ 13.3, T ■ 2.0. 
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Fig. 13. Coapariaor of obacrvad (OHAC) and thaoraclcal tlae hiatorlta for the 
aagnltude of the aagnetlc field for Pionaar Venus orbit 6 baaed on inbound and 
outbound interplanetary conditions using gasdynaalc solution for K. “ 13.3, V ■ 2.0. 
(a) Magnetic field nagnltude. (b) Magnetic field coaponents. 
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Fig. 14. Plonaer Venua orbit 3 trajtccorlaa and obaarvatlonal bow ahock croaainga as 
viewed In solar wind coordinates based on Inbound and outbound Interplanetary solar 
wind directions; also, various bow shock shapes for different interplanetary solar 
wind conditions. 


netlc field. With the known (Bx • By , Bz ) 
components the corresponding components (Bj^, 
®v. • rotated system are given by 


fB' 1 1 

1 i 1 

fi 

0 

” ! 

'°y*i ' ' 

0 

1 

cos9 

sinS ' 

1 

|b; I 1 

r 

-sln6 

cos6 1 
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1 i’b„ 


(52) 


r^-j 


In this reference frame the parallel, perpendic- 
ular, and normal Interplanetary components are 
identified as 


B' 

X 


B' 

V 


B' 

z 


(53) 


a ■ tan 
n 


-1 


B 


/(B_ )2 + (B. )•' 


■ tan 


-1 


B' 

z 


.^(B* ) + (B' )^ 


(55) 


The magnetic angle lit associated with the Incident 
perpendicular component and the unit magnetic 
field ratios (|b|/B_)„, (|||/B^)^, and (B/B^)„ In 
the rotated system are next determined by bilinear 
interpolation through the flow field grid. Then 
the magnetic field components (B^, By, Bj) In the 
rotated system are calculated from 


B' 

X 


+ cos* 


cos* 


slna' 

P 


P 

I B I 

i“!. 


i - 

|i^ 


(56) 


The magnetic field angles a' and a* in the rotated 
system are given by ** " 


B’ - coso' 

y n 


sin* 



[b 1 


’ B* 

1 

1 § 


tan-' 

r 

00 

• • 

1 

- tan-' 

y. 

B' 

L *-J 

(54) • ‘inOp • 



B 

b” 


( 57 ) 
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Fig. IS. CoBparlson of obaarvcd and thaoretlcal time hlstorlei of lonoaheach plaana 
properties for Pioneer Venus orbit 3 based on Inbound and outbound Interplanetary 
solar wind rondltlons. 



B' • slna’ • • B (58) 

The magnetic field cooponents In the solar wind 
(X, T, Z) system are Chen determined from the 
rotational transformation 


with B^ either parallel or perpendicular to y^. 
These'and similar decermlr.atlons for the following 
35 other sets of conditions have been Included In 
a catalog of results available from Stahara et al. 
[1980]: 

- (2.0, 3.0, 5,0, 8.0, 12.0, 25.0) 


■b 1 fi 0 0 I |b;1 

By • 0 ccs6 -sine B^ ! (59) 

B^ 0 sine cose B^ 

and finally In the sun-planet system from 




-cosOcos^ sln.'i -cos.lsln# [By 

P P * 

slnOcost -cosO slnflslnd B„ (60) 

P P ‘ 


Rasults 

Using the computational procedures described 
above, a large number of solutions have been 
calculated for a wide variety of conditions repre- 
sentative of those that might be experienced at 
Venus. A swple of these results is illustrated 
In Figure 9 for ■ 8, 7 • 5/3, and H/Rj, ■ 0.10 


H/R - (0.01, 0.10, 0.25) 
o 

H/R - (0.10, 0.20, 0.25) 
o 

Also Included in that reference are a few results 
for conditions well outside the range likely to be 
enco«Bitered at Venus but Chat are useful beesuse 
they Illustrate the extended capability of Che 
present procedures. 

Verification of the new procedures has been 
done In a variety of ways In addition to compari- 
son with previous results for sets of conditions 
for which both sets of results could be determined. 
These Include consideration of a variety of 
special test cases In which the location In the 
flow field and the Incident Interplanetary magnetic 
field orientation were systemaclc^ly changed so 
as CO produce both sywetrlc and antisymmetric 
changes In the resultant lonoaheach magnetic 
field, as well as to reverse the roles of the 
perpendicular and normal components. All of these 
various permutations of the magnetic field calcu- 
lations were successfully verified. 

The final and ultimate check of Che theoretical 


i 
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Fig 16. Conparlton of observed (OKAO) and theoretical tine histories for the magnetic 
field for Pioneer Venus orbit 3 based on inbound and outbound Interplanetary solar wind 
conditions using gssdynanlc solutions ■ 7. 33, F ■ 2.0 for Inbound and “ 5.96, 

7 • 2.0 for outbound calculations, (a) Magnetic field aagnltudes. (b) Magnetic field 
coaponents. 


nodel and nvaerlcal procedures lies In the coo- 
parlson of the predicted results with data 
actually aeasured by a spscecrsft. Attention Is 
now turned, accordingly, to the results of pre- 
1 Ixlnary cosparlscns with data obtained fron 
orbits 3 and 6 of the Ploneer-Venus orblter. 

These orbits were selected In conjunction with the 
experlse Iters on the basis of completeness of data 
and apparent steadiness of conditions during the 
time of Interest for the study. 

To perform calculations suitable for comparison 
with data from a particular orbit. It Is necessary 
to specify appropriate conditions corresponding to 
those actually measured In both the Ionosphere and 
In the oncoming Interplanetary solar wind plasma. 

On the basis of data from the orblter retarding 
potential analyzer (OPKA) Indicating the Iono- 
spheric scale height to be approximately 200 km 
for orbits 3_and 6, a value of 0.03 was selected 
for H/Rq or H/Kg. For such small values the two 
Ionospheric pressure models Indicated by (20) and 
(21) yield essentially the same obstacle shape, 
as can be seen from Figure 2. With regard to on- 
coming Interplanetary conditions It Is required 
that values for the solar wind bulk velocity y_, 
density o^, temperature T_, and magnetic field 
be specified. The first three are provided by 
the orblter plasma analyzer (OPA), and the 
magnetic field Is given by the orblter flux gate 
magnetometer (OHAC) . OPA provides either Ion or 
electron data, but not both simultaneously. For 
orbits 3 and 6, Ion measurements were available 
and have been employed. Data on the oncoming 
direction of the solar wind, as given In terms of 
the angles (Q, dp), define the coordinate rotations 
required to align the gasdynamlc calculation In the 


free-stream solar wind direction; while Information 
of solar wind speed, density, and temperature serve 
to define the free-stream gasdynamlc Mach number 
M^. 

Although it Is known the temperatures of the 
Ions and electrons may be substantially different 
In the solar wind plasma (a ratio - 4 Is 

often representative), they are assimMd to be equal 
(T«/Ti “ 1) In the derivation of the single-fluid 
magnetohydrodynamlc equations upon which the 
present calculations are based. Now that the 
theory has been put Into a more readily usable 
form and more extensive data on the temperature 
of both zlectrons and Ions are becoming available 
for comparative studies, steps should be taken to 
Improve understanding of the relation between 
single- and multiple-component plasma theories for 
space plasma flows, and the calculatlve procedures 
refined as appropriate. 

With this Information the detailed gasdynamlc 
and unit magnetic field calculatlona In the lono- 
sheath region can be carried out. To provide an 
Idea of the detail obtained by the present com- 
putstlonal procedures. Figure 10 displays the flow 
field grid for one of the gasdynamlc solutions 
used In the comparisons. Values for each of the 
plssma propertlea y/v^, o/o„ T/T., (5 /B_)m, 

(B/B )x and (B/B^)_ are determined at each Inter- 
section of the grid lines, including the bow 
shock, stagnation streHllna, and lonopause 
boundary. The final output of the calculation 
consists of these quantities and the coordinates 
of the bow shock and lonopause, expressed either 
In tabulsr form or as plotted contours, and also 
as time histories along a specified trajectory. 

For comparisons with spacecraft data the time- 
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Fig. 16. (continued) 


history predictions ere the uost convenient. For 
the theoretical predictions, two sets of results 
are usually generated on the basis of the last/ 
first Interplanetary solar wind properties 
■easured before./af ter bow shock inbound /outbound 
crossing. 

Figure 11 displays sooe overall flow field 
results for orbit 6. Indicated are bow shock 
locations for the three coabinatlons of and 
plasma specific heat ratio 1 , l.e., (M^, >) - 
(13.3, 5/3), (13.3^ 2), (3.0, 5/3) for flow about 
an ionopauae with H/Rq • 0.03. Also Indicated are 
two sets of points (solid lines with circles, and 
dashed lines with squares) representing the space- 
ertft trajectory as viewed in two solar wind 
oriented coordinate systems. The trajectory 
indicated by the solid lines and circles is that 
based on the last measured direction ( 0 , gp) ■ 
(6.5*, -1.4*) of the interplanetary solar wind 
Just prior to crossing the bow shock on the in- 
bound leg, while the dashed lines and squares 
denote the trajectory based on the first measured 
direction (f,, Sp) - (4.9*, 7.6*) of the solar wind 
iasaediately after crossing the bow shock on the 
outbound leg. These results Illustrate the 
extremely large dependance of apatial position of 


a trajectory point, as viewed in solar vlnd coor- 
dinates, on solar wind direction. For the partic- 
ular angles Indicated the shift in X coordinate 
of a trajectory polr.. can be as great as a quarter 
of the Venusian planetary radluj and obviously 
results in substantial differences in predicted 
flow and magnetic field properties. Although 
Spreiter and Rlzcl 11972] noted the Importance of 
such an effect in the interpretation of the data 
from Mariner 4 in its flyby of Mars, the effects 
of angular shifts in the solar wind direction are 
usually Ignored in most discussions of the inter- 
action. The present results, however, indicate 
that this purely geometrical affect can be very 
slgnlf leant, even for directional shifts of less 
than 5*, and must be accounted for in any 
quantitative comparison of theoretical and 
observational results. 

The selection of values for M^ and Y in 
rigure 11 represents an attempt to resolve the 
uncertainty in the appropriate values for these 
parameters. Because only solar wind ion temper- 
atures from the OPA were available for orbit 6, 
the initial calculations of were based on the 
assumptions thst T, - T^, which loads to • 13.3 
and Y - 5/3, as appropriate for a gas with parti- 
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Fig. 17. Comparison of observed (OMAG) and theoretical time histories of the magnetic 
field for Pioneer Venus orbit 3 based on Inbound solar wind Interplanetary conditions 
using a gasdynamlc solution for • 3.0, Y • 5/3. (a) Magnetic field magnitude. 

(b) Magnetic field components. 


cles possessing 3 degrees of freedom for their 
molecular motion. The resulting location for the 
shock wave Is not in good agreesient with the 
observations. To Investigate whether changing 
the value for Y to 2, which corresponds to a gas 
with 2 degrees of freedom ae might be supposed to 
occur If the magnetic field sufficiently aligns 
the molecular motion, the calculations were 
repeated with M^ “ 13.3 and Y > 2. The results. 
Indicated by the dashed line, are In Improved, 
but still not completely satisfactory, agreement 
with the observed shock locations. Finally, If 
It Is assumed that the oncoming Interplanetary 
electron temperature Is not equal to the Ion tem- 
perature but Is substantially higher, we are led 
to consider low values of the order of 3-5 for M . 
Results for M^ ■ 3.0, Y ■ 5/3 are Indicated In ” 
Figure 11 by the solid line. It may be seen that 
results calculated for these valuas and for the 
Inbound solar wind direction are In very good 
agreement with the theoretical results. These 
comparisons make evident both the possibilities 
and the uncertainties In the use of the theory at 
the present level of understanding of the parameter 
selection procesa. 

Figure 12 displays the tlsM history comparisons 
of the predicted bulk plasma density, speed, and 
temperature with OPA data. The solid lines with 
circles correspond to results based on Inbound 
interplanetary conditions, while the dashed lines 
with squares correspond to outbound conditions, 
xhile the one data point in the lonosheath for 
each quantity Is In general agreement with the 
theoretical calculation, the lack of more detailed 
plasma data prevents a definitive conclusion. The 
OPA requires approximately 9 min to acqul-e suffi- 


cient data to enable determination of the bulk 
plasma quantities. This time Interval usually 
presents no problem when the spacecraft Is In the 
Interplanetary solar wind, but the large resolu- 
tion time effectively averages the plasma quan- 
tities In the lonosheath over such a large spatial 
range that only overall comparisons of the bulk 
plasma properties are possible with this 
Instrument. 

The situation Is quite different for the 
magnetic field, since the OMAG provides essentially 
Instantaneous measurements. Figure 13 displays a 
comparison of the data from orbit 6 with two sets 
of theoretical results based on the Inbound and 
outbound Interplanetary conditions. In these com- 
parisons both sets of predictions exhibit good 
agreement with the data for the appropriate part 
of the trajectory. Thus, on the Inbound leg the 
predictions based on the Inbound Interplanetary 
conditions are In good agreement with the data, 
while the outbound-condition predictions are 
notably inferior, particularly with regard to 
shock crossing; and conversely. Since steadiness 
of the Interplanetary conditions was a criterion 
used In the selection of orbit 6 for these com- 
parisons, the results of Figure 13 serve to 
Illustrate the need for simultaneous measurements 
of Interplanetary conditions and planetary flow 
field conditions in more definitive evaluations 
of the theory. 

Results for corresponding comnarlsons for 
orbit 3 are given In Figures ] ', and 16. In 

Figure 14, bow shock location ^resented for 

five different coaiblnatlons oi , .<nd Y. The Mach 
numbers • 7.38 and 5.96 correspond to the In- 
bound and outbound Interplanetary conditions. 
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Fig. 17. (continued) 


respectively, for c^, and T as measured for 

the ions by the OPA; the two values for ' " 5/3, 

2 are indicative o/ present uncertainty about the 
effective ratio of specific heats for the solar 
wind plasma. The bow shock location for M,, * 3 
and Y - 5/3 given previously In Figure 11 for 
orbit 6 Is also Included, because once again it is 
closer to the observed shock crossings. Also 
provided In Figure 14 are the orbital trajectories 
as viewed in solar wind coordinates for the In- 
bound (il, ♦„) ■ (3.3*, 0.15*) and outbound (SI, tp) 
- (3.7*, 4.9*) solar wind directions. 

The comparisons for the bulk plasma properties 
for orbit 3 shown in Figure 15 are urked again by 
the paucity of observational data. What informa- 
tion there Is suggests agreement for the plasma 
speed and density but a notable discrepancy for 
the temperature. This is tho.^.it to be indicative 
of the manner in which the bulk properties from 
the theoretical model are interpreted in relation 
to the observational meaauremevts; l.e., the theo- 
retical values correspond to those for a single- 
component plasma in which it is considered that 
the electrons and ions have equal temperatures, 
whereas the measurements are for the ions only in 
a multicomponent plasma in which the ion and 


electron temperatures are usually quite different. 

A more penetrating analysis of the meaning of the 
temperature in the theory and its relation to that 
deduced by the axperlmenters from their data 
appears to be a necessary and Important subject 
for future study. 

Results for the magnetic field comparisons are 
displayed in Figure 16, which provides time 
histories of the magnitude and individual com- 
ponents based on both Inbound and outbound inter- 
planetary conditions. Although the shock crossings 
are somewhat in disagreement, since the calculated 
results are for M,, ■ 7.56, 6.96 and Y - 2, the 
general trends displayed by the calculated and 
observed results are notably similar, particularly 
if account is taken of the drift with time along 
the trajectory of the transition from the calcu- 
lations based on inoound conditions to ones based 
on the outbound conditions. 

To demonstrate the improvement obtained in 
magnetic-field results when a gasdynamic flow 
field solution is employed which more closely 
agrees with the obaervational bow shock location, 
the magnetic field time histories from orbit 3 are 
cooq>ared in Figure 17 using gasdynamic results for 

■ 3.0, Y • 5/3, and the directions (H, «p) • 
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Fig. 18. Illuscraclon of effect on bow ab./ck location of variability of oncoBlng 
solar wind direction; M • 8.0, > 5/3 flow past a sphere /cylinder lonopause 

obstacle. 


(3.3°, 0.15°) measured at the Inbound shock cross- 
ing. As can be seen, there Is a narked Improvement 
In the agreement near the bow shock and quite good 
agreement throughout the lonosheath for borh the 
magnitude and the Individual magnetic field com- 
ponents. In addition to illustrating again that 
further study needs to be made of the parameter 
selection process required to relate the calculated 
and observational results optimally, these compari- 
sons emphasize once again the need to account Ir. 
the calculations for the actual direction of the 
Interplanetary solar wind flow. 

To Illustrate the significance of this point. 
Figure 18 has been prepared to show that even a 
modest change of ±6* In the direction of the 
Incident solar wind can have substant^::! elfects 
on the location of the lonopause »':a bow wave, as 
viewed In a coordinate system 'signed with the 
Sun-Venus direction which 1'. typical of most 
presentations of observational data. It would be 
of considerable Interest to determine how much of 
the scatter of data for the location of the bow 
wave may be attributable to this simple factor. 

A further point of growing concern Is assocl ■ 
ated with effects of significant changes down- 
stream of the terminator of the lonopause shape 
from that calculated herein. More detailed :cm- 
parlsons will have to be reserved for future 
studies, fut the resw'^ts of Figure 19 for the 
theoretical characte. stlc lines emanating from 
the Intersection of the lonopause and the termi- 
nator define the region of Influence of a small 
Inward or outward caper of the lonopause profile 
downstream of the terminator. For the entire 
region displayed extending dosmacream to more chan 
16 Venus radii, these results Indicate that there 
Is no accompanying change in the location of the 
bow wave. It may also be noted chat Che point 
designated by nu^er 2 In the Mariner 5 data from 
Venus (see Spreiter and Alksne 11970] for an 
earlier comparison with theory), where the lono- 
sheath plasma changes from comparatively quiet 


CO disturbed. Is romarkably close to the charac- 
teristic line from the cerm.'nacor. Such behavior 
could result from unsteady fluctuations of the 
lonopause surface beginning approximately at the 
terminator. The data from Pioneer Venui should be 
examlneo to determine whether this Is a general 
property of the plasma In Che vicinity of these 
characteristic lines or perhaps Just the result 
of some passing transient occurring at the time 
of the Mariner S encounter. 


Concluding Remarks 

The present application of advanced computa- 
tional procedures to the modeling of solar wind 
flow past Venus was undertaken to Improve the 
accuracy and utility of the theoretical predic- 
tions. Starting with the steady, dlsslpaclonlkss, 
magnetohydrodynamic model for axlsyMecric, super- 
sonic, super-Alfvdnlc solar wind past a nonmagnetic 
planet with a shielding ionosphere, a nuzd>et of 
Important theoretical extensions have been devel- 
oped and Implemented. These Include the capability 
for creating lower Interplanetary gaadynamic Mach 
numbers M_ down to about 2 and a wider variety of 
lonopause shapes Including a new family of ahapes 
which Includes the effect of gravitational varia- 
tion in the scale height. Additionally, the 
capability for determining the plasma gaadynamic 
and magnetic fle'd properties along any arbitrary 
trajectory, accounting for an arbitrary oncoming 
direction of Che solar wind, has been established. 
All of these dcvelopmencc havo been Incorporated 
Into a current computational model to enable 
detailed calculations of Che solar wind Inter- 
action with Vanus. The model has been compre- 
hensively described by Stahara et al. [1980]. 

Also Included In that account Ic an axtensive 
catalog of results for a wide range of flow con- 
ditions and lonopause shapes representative of 
Chose chat might be anticipated at Venus or at 
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Fig. 19. Illuatratlon of effect of various wake shapes on bow shock location and 
characteristics from terminator; Y - 5/3, ■ 3 and 8. 


some other nonmagnetic planet with a shielding 
ionosphere . 

Finally, theoretical results determined from 
the model have been compared with data from two 
of the early orbits of the Pioneer Venus orbiter 
These comparisons have indicated the Importance, 
heretofore largely neglected, of the directional 
variability of the oncoming solar wind. These 
results, taken In toto, serve to verify the basic 
theoretical model and the effectiveness of the 
numerical solution, although It is evident that 
further study should be applied to certain 
questions relating to Interpretation of plasma 
quantities like the temperature. The results 
demonstrate also the value of the present theo- 
retical procedures as a research tool capable of 
routinely providing — at small computational cost 
and in a form directly compat'ible with observa- 
tions — details of the solar wind/planetary 
atmosphere Interaction process not previously 
attainable. 

With regard to future uses as well as improve- 
ments of t> e cresent model, the obvlrjs need for 
a detailed st idy Involving comparison between 
theory and observat.ions for a large number of 
orbits of the Pioneer Verjs orbiter is clear. 

Or. the basis of the preliminary comperlsons for 
'■ oits 3 and 8 the calculated sLaguetic field 
apoears to be remarkably accurate for relatively 


qulettlme conditions. Similar comparisons of the 
plasma properties Indicate a need for an Improved 
Interpretation of the results from the single-fluid 
theory in terms of multiconponent measurements. 
Questions regarding the possible suppression by 
the Interplanetary siagnetic field of the number 
of degrees of freedom of the plasma require 
further study and could be clarified through 
systematic conparisona with data. Additionally, 
observations from the Pioneer Venus orbiter of 
the nightside Ionosphere of Venus have revealed 
a more complex and dynamic structure than sus- 
pected. These observations point. In particular, 
toward the need for Improvement of Che simple 
model lised In the present method for Che deter- 
mination of the Ionosphere boundary downstream of 
the terminator. This improved determination would 
involve an Iterative procedure In which a balance 
of Che sum of the solar wind gasdynamic and 
magnetic pressures along the ionopause surface 
would be maintained against Che Ionospheric pres- 
sure, not necessarily assumed to be spherically 
symmetric The present method, which balances 
the extended Newtonian gas plus magnetic pressure 
distribution against the Ionospheric pteesure, 
might represent the first step of such ar. Itera- 
tion; but it Is also possible that dynamic effects 
associated with Intermittently passing volumes of 
magnetized plasma tom from more upstream regions 
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of cha lonopauaa nltht (raacly Impair tha ganarai 
utility of any atatlonary modal for tha lonopauaa 
downatraam of tha tarmlnator. Thla la anothar 
Important topic awaiting furthar atudy. 
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